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Abstract 

The leading nonperturbative effects to the fine and hyperfine splitting were 
calculated some time ago. Recently, they have been used in order to obtain 
realistic numerical results for the lower levels in bottomonium systems. We point 
out that a contribution of the same order 0{h.Q CD /nrt'a^) has been overlooked. 
We calculate it in this paper. 
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1 Introduction 



The leading non-perturbative (NP) corrections to the fine and hyperfine splitting 
were calculated a few years ago with the theory developed by Voloshin and Leutwyler 

for heavy quarkonia. Leutwyler already computed the hyperfine splitting for 
n — 1, I — in ref. [fjj. However, a contribution of the same order 0(AQ CD /m 3 a 2 s ) was 
overlooked. Later on, a number of authors (Voloshin, Di Giacomo et al. and Dosch 
and collaborators) calculated the fine and hyperfine splittings for different quantum 
numbers (n — 2, 1; I — 0, 1) in refs. || where the above overlooked contribution was 
pointed out though finally neglected. 

Recently, new efforts have been made to obtain realistic numerical results from 
rigorous QCD in refs. [f|, |[ where a coherent picture of bottomonium with n = 1, 2 is 
provided. Relativistic, radiative and the leading NP corrections were put together for 
the first time. Unfortunately, the above mention contribution was also omitted. 

The aim of this paper is twofold: (i) to complete the calculation up to order 
0{K , Q CD /m?a 2 s ) and (ii) to control the error made by neglecting this contribution. 
The last point is important if we eventually would like to improve these results by tak- 
ing into account further orders in the perturbative and NP expansion. We also provide 
results for a wide range of quantum numbers. 

Let us briefly discuss how this contribution arises. Consider the Breit-Fermi inter- 
action in the octet potential. 

yCoul _^ yCoul + yBF ^ ^ 

then the octet propagator changes into 

1 (i2) 



E — Hg E — — Vg* F E — H 8 E — 8 E — Hg 

where H 8 is the octet Coulomb Hamiltonian. This leads to the following correction to 
the energy 

OiV c £j n — li 8 £j n — H 8 

E n is the Coulomb singlet bound state energy, n is the principal quantum number. I, 
s and j are the angular momentum, spin and total angular momentum, respectively. 
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Let us comment upon the relative size of this contribution. Notice that it is 1/N C 
suppressed (in fact jv2 1 _ 1 ) by the ratio of octet 0(^-) to singlet 0( i ^" 1 ) coupling 
constant and also it is suppressed by large energy differences between octet and singlet 
states (recall that the octet potential is repulsive). From this we can conclude that 
matrix elements with a larger number of octet propagators are numerically suppressed. 

Notice also that V^ F only affects the leading NP results for the fine and hyperfine 
splitting. For the standard NP results found by Leutwyler and Voloshin in ref. |IJ |2| 
( |1.3|) gives rise to a subleading contribution which will be neglected in the following. 
Therefore, we only take into account the piece of V^ F which contributes to the fine 
and hyperfine splitting 



vr = v? + vr = v 8 ls + vf + v» f (i.4) 



where 



S u (r) = 2(3S 1 -fS 2 -r-S 1 -S 2 ), (1.6) 



and 



5E 8 = 5E F + 5E^ F = 5E^ S + 5E F + 5E^ F . (1.8) 

We shall give analytical formulas for the above energies when available. 

Let us define the physical fine and hyperfine splitting for n — 2, I = 1. We use 
spectroscopic notation 2 2s+1 Pj. The generalization to arbitrary quantum numbers is 
straightforward. The fine splittings read 

A 21 = M(2 3 P 2 ) - M(2 z P l ) , Aio = M(2 3 P!) - M(2 3 P ) (1.9) 

while the hyperfine splitting reads 

A HF = M(2 3 P) - M(2 1 P 1 ) (1.10) 

where 

M(2 3 P) = i(5M(2 3 P 2 ) + 3M(2 3 Pi) + M(2 3 P )) (1.11) 
9 

Typically, 5E HF only contributes to A HF and SE F only does to A 2 i and A 10 . Certainly, 
in our case 5E% F only contributes to Ahf and A 2 i, Aio only get contributions from 
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8E[ but rather peculiarly, as we will see below, A HF receives contributions from 5E$ , 
that is, 5 E[ ^ where 

6E?= (1.12) 
{(2/ + 3)^(j = J + 1) + (2 J + l)^(j = + (2i " 1)^0' = * - !)} 

for a general /. This is quite unusual and did not happen with the NP corrections 
calculated previously in the literature where 5E F = 0. 



2 Fine Splitting 

For the fine splitting we obtain 

We have split the matrix elements in radial and angular integrals. 

Ill 

F(n,l,h,l 2 ) = (R n i\r _ {h) ^ _ H (i 2 ) r \ R m) ( 2 - 3 ) 

ti 8 ^8 

and 

G LS (j,l;h,l 2 )= £ £ (ai,« = l)|r i /M^Wiia« 1 « = l)) (2-4) 

jl = |h-l|j2 = K2-l| 

G T (j,/;/i,/ 2 )= £ £ <0\^« = l)l^5 12 (0W i |0',/,a = l)> (2.5) 

Jl = Kl-l|j2 = |«2-l| 

where R n i(r) is the radial Coulomb wave function, 

\(n,j,l,s)) = \R nl )\(j,l,s)), (2.6) 

1 10= - * ,J Q, (2-7) 



and 7j ; ; is the identity in the subspace with total angular momentum j, orbital mo- 
mentum I and s = 1 (otherwise the matrix element is 0). 
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Eq. ( |2.3| ) can be computed using the techniques shown in ref. || but we have 
not succeeded in finding a close analytical expression for F(n, Z; Zi, I2) although we did 
succeed for the angular momentum functions. They read 



G Lb (j,l;h,l 2 ) = (-iy +h+1 5 h>h (2l + l)C(l,l,h;00y 



G T (j, Z; l u h) = (2/ + 1)C(Z, 1, Z l5 00)C(Z, 1, Z 2 ; 00) 



ii+i 



x 



h ji 1 1 f Z2 ji 1 



E (% + !) • , , (2-9) 



Jl = |/l — 1| 



J z 1 I j / 1 



x2 (5 W2 - 3C0i, 1, Zi; 00)C(ji, 1, Z 2 ; 00)) . 

We display the explicit expressions in two tables. 

Finally, the final expressions for the energy corrections read 



5E L S S = 5.,^ * { % G \ fls[n,l,j] , (2.10) 



a s ii{a s G 2 



5E T S = 8 S ^ ^ G X ft[n,l,j] . (2.11) 



s J 

a s was defined in ref. M. 

We write ( 2.10j) and ( |2.11| ) in this way in order to make comparison with || simpler. 



We give some numbers for lower values of n. 



(fls + ft) [1, 0, 1] = 9868 , (fls + ft) [2, 0,1]= 236464 , 

\J 8128125 U 19780605 

p , rnin1 15475732 . ri „ Nro 1 . 8134524806 

(/Zs + ft) 2, 1,2 = , (fls + ft) 3, 1,2 = 

u ^ n ' ' J 415392705' U ^ n ' ' J 53682451515 

Mrn , -452188 ,„ „, rni „ n 6745478 

{fls + /t)[2, l, 1] = — — — , (fls + ft)[3, 1, 1] 



19780605 ' w 59449005 ' 

1 -1235984 ,„ „ Nr n 83789896 
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In principle our contributions for the fine splitting are quite small in comparison 
with the NP and radiative corrections coming from the wave functions. Nevertheless, 
the authors of ref. || managed to isolate the latter in an unknown factor which is 
determined from data. Hence, the remaining perturbative and NP effects are kept 
under control. For the fine splitting they obtain 

A 10 = - (1 + 8 rad ) A 21 - 5 NP . (2. 13) 

They take A21 from data and 8 ra d and Snp are respectively the remaining radiative 
and NP corrections (see (3.2) in the second paper of || for details). While, the new 
contribution reads (to be added to 5np in ( [2.13| )) 



«. / x - 5 r^F/ x 9 . ^f, , a s n(a s G 2 ) 2548892 

SMnew) = - SEib = i)+l5Eib = iysES(J = o) = -^ m ; (Cf 4, 415392705 

(2.14) 

being around 1% smaller. 



3 Hyperfine Splitting 

In this section we work out the analytical formula for the hyperfine splitting. It 
receives contributions from both SE F and 8E% F . Let us calculate SEg F . The angular 
integral is almost trivial. The delta function only survives for angular momentum 
which after composing with r implies that only I = 1 states contribute. It is quite 
remarkable that only for I = 1 states we obtain a non-zero contribution. After these 
manipulations we find that SE$ F becomes 

8E£ F = s{s + m JL as J^ sG2) \{R nl \r l -n-x\r = 0)| 2 . (3.1) 

8 v J ' 2N c 9m 2 6N C IN 1 E n - H^ 1=0) 

The radial integral can be carried out by following the formulas given in 0. We obtain 

6E* F = 6 ltl s(s + l)a 8 n } asG *h f[n] , (3.2) 

n 5 (n 2 -l)r[9n/8-2] 2 
hf[H] '= 2* r[9n/8 + 3] 2 (3 ' 3) 
We again write ( |3.2| ) in this way in order to make comparison with || simpler. We 
give some numbers for lower values of n (n > 2). 
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[2] = ^5536 J 07772, „ 



32967675 ' 296815671075 ' 

hftt]- 2048 hm= 2097152000 

J 135270135' ■*{ \ 324918632936187' 
Let us know discuss our contribution to Ahf- If / = it turns out to be quite small 

but for I the leading perturbative order is 0. Therefore, the next perturbative order 
is needed and in principle the nonperturbative effects are going to be more important. 
Thus, in the second paper of a careful analysis of the hyperfine splitting was done 
for n = 2, I = 1. It was obtained 

5 p 21 976 7r(a s G 2 ) 

Ahf = 24 ( T " T )CFaA * + T0E3 m%C F a s y " (3 " 5) 
While our contribution reads (to be added to (|3.5| )) 

A ( \ JW^^hf a s n(a s G 2 ) 160277456 

A HF (new) = 6E[ + 5E 8 = ^ m3( ^ )2 18692671725 ( 3 - 6 ) 

which turns out to be around 1% smaller. 



4 Discussion 

From the phenomenological point of view our results are going to be important 
only in the event that a very high precision measurement is done. Nevertheless, they 
are conceptually important since they take into account the perturbative corrections 
to the octet Coulomb potential. We have also seen that the hyperfine splitting gets 
corrections from the fine octet potential which is something somewhat unusual. 

Our results complete the calculation of all the contributions of order 0{h.q CD /w?a 2 s ). 
We have also seen the error in neglecting this contribution. This step is unavoidable 
for an eventual improvement of these results by calculating higher orders in the per- 
turbative and NP expansion. 
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J =1-1 


3 = 1 


J =1 + 1 


h = l-l 


l-l 2 
1+21 


l-l 

1+21 


(-1+1)1 
1+21 


h = I 











h = l + l 


(1+0 (2+0 
1+21 


2+1 
1+21 


1(2+1) 
1+21 



Table 1: We display here G LS . For I — j must be 1. The remaining matrix elements 
are zero. 



j = 1-1 



l 1 = l- l,l 2 = l- l 



-6-111-261 



2 -Z 3 +10 1 4 -12 1 5 -8 1 6 ) ^ 



I (—3+2 Q ( — 1+2 Q (1+2 1) 



2 (-1+0 (2+3/+8/ 2 +4/ 3 
I (-1+21) (1+2 1)' 2 



h = l-l,l 2 = l + l 
-6(1+0 



(1+2 iy 



l 1 = l + lj 2 = l + l 

2 (1-2 Q (1+Q (2+Q 

2 (-6-5Z+7 1 2 +12 1 3 +4 1 4 ) 
(1+2 I) 2 (3+5 1+2 1 2 ) 



3 = 1 



d+2 iy 



2 (-6+14/+37« 2 +« 3 -50 1 4 -36 I s -8 / e 
(1+Q (1+2 Q 3 (3+2 Q (5+2 Q 



3 = 1 + 1 



2 (l-Q Z (3+2 Q 
(1+2Q 3 



-6Z 
(1+2 Q 3 



Table 2: We display here G T . For I = j must be 1. The remaining matrix elements 
are either zero or they can be deduced by exchanging l\ ++ l 2 . The asterisk indicates 
that the result is only valid for / > 2 otherwise it is 0). 



